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A study of the conductance noise in a two-dimensional electron system (2DES) in Si at low 
temperatures (T) reveals the onset of large, non- Gaussian noise after cooling from an equilibrium 
state at a high T with a fixed carrier density Ug. This behavior, which signifies the falling out of 
equilibrium of the 2DES as T ^ 0, is observed for Us < Ug [ug - glass transition density). A 
protocol where density is changed by a small value Aris at low T produces the same results for the 
noise power spectra. However, a detailed analysis of the non- Gaussian probability density functions 
(PDFs) of the fluctuations reveals that An^ has a qualitatively different and more dramatic effect 
than AT, suggesting that Aris induces strong changes in the free energy landscape of the system 
as a result of Coulomb interactions. The results from a third, waiting-time {tw) protocol, where ris 
is changed temporarily during tw by a large amount, demonstrate that non-Gaussian PDFs exhibit 
history dependence and an evolution towards a Gaussian distribution as the system ages and slowly 
approaches equilibrium. By calculating the power spectra and higher-order statistics for the noise 
measured over a wide range of the applied voltage bias, it is established that the non- Gaussian noise 
is observed in the regime of Ohmic or linear response, i.e. that it is not caused by the applied bias. 

PACS numbers: 71.55.Jv, 72.70.+m, 71.30.+h 



I. INTRODUCTION 

Many novel materials find themselves close to the 
metal-insulator transition (MIT), where strong electron- 
electron interactions and disorder are believed to give 
rise to a variety of complex phenomenal Near the MIT, 
these two effects are usually comparable in magnitude, 
so that their competition is expected to lead to glassy 
behavior of electrons, in analogy with other frustrated 
systems j^'*^ A common denominator for all glasses is the 
existence of a complex or "rugged" free energy landscape 
(PEL), consisting of a large number of metastable states, 
separated by barriers of different heights. This results 
in phenomena such as slow, nonexponential relaxations, 
divergence of the equilibration time, and breaking of er- 
godicity, i.e. the inability of the system to equilibrate 
on experimental time scales. Therefore, such out-of- 
equilibrium systems also exhibit aging effects, where 
the response to an external excitation (i.e. relaxation) 
depends on the system history in addition to the time t. 
A detailed analysis of temporal fluctuations (noise) of the 
relevant observables yields complementary information 
on configurational rearrangements or transitions between 
metastable states. Non-Gaussian distributions of various 
observables in glassy systems have been reported,^ re- 
flecting the presence of large, collective rearrangements. 

Most experimental studies of charge or Coulomb 
glasses have focused on situations where electrons are 
strongly localized due to disorder, far from the MIT.^ In 
recent years, however, a two-dimensional electron system 
(2DES) in Si metal-oxide-semiconductor field-effect tran- 
sistors (MOSFETs) has emerged as an excellent model 
system for studying glassy or out-of-equilibrium charge 
dynamics near the MIT.^ In particular, studies of both 
relaxations and noise on very disordered samples have 



established that the 2DES in Si exhibits all the main 
manifestations of glassiness: slow, correlated dynamics 
(non-Gaussian noise), ^^^^ nonexponential relaxations^ 
diverging equilibration time (as temperature T ^ 0)^^ 
aging and memory. — ~— Glassiness is observed for all car- 
rier densities Ug < Ug, where Ug > ric {ric - the critical 
density for the MIT), thus giving rise to an intermediate 
glassy, "bad metal" phase {kpl < 1; - Fermi wave vec- 
tor, / - mean free path). In a 2DES with a relatively low 
amount of disorder, noise studies show that the interme- 
diate phase becomes vanishingly small, i.e. Ug > Uc-^^'^^^ 
These observations are consistent with predictions of the 
theoretical models that describe the MIT as a Mott tran- 
sition with disorder. The results presented here were 
obtained on samples with a relatively large amount of 
disorder, where the intermediate {ric < Us < Ug) glassy 
phase is more pronounced, as mentioned above. We note 
that aging studies have found^^^^^ that an abrupt change 
in the nature of the glassy phase occurs at nc, i.e. at the 
2D MIT, signifying the intrinsic glassiness of the 2DES. 

Glassy dynamics in a 2DES in Si has been studied us- 
ing several different experimental protocols, all of which 
involved rapid changes of Ug ("Arig protocol"). For ex- 
ample, a large, rapid (within 1 s) change of at a low 
enough T results in a nonexponential relaxation of con- 
ductivity ait) towards a new equilibrium value. Small 
changes of Ug at low T do not produce any observable 
relaxations,^^ but have led to the observations of non- 
Gaussian fluctuations in (j{t)S'^-^ On the other hand, 
in other types of glasses (e.g. structural, spin), studies of 
aging processes and other glassy properties are typically 
carried out following a thermal quench from an equilib- 
rium state at a high T to a nonequilibrium state at a low 
T. The analysis of the relaxations following such a ther- 
mal quench, however, is greatly complicated by their non- 
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trivial dependence on the cooling time and, hence, a fairly 
fast cooling process is essential (see, e.g., Ref. 17'). It is 
of fundamental interest to establish whether the charge 
glass realized in a 2DES in Si also falls out of equilibrium 
as T is reduced, in common to other types of glasses. 

Here we investigate the charge dynamics in a 2DES 
in Si at low n^, near the MIT, after cooling from an 
equilibrium state at a high T with a fixed ( "AT pro- 
tocol"). Since it is not possible to achieve rapid cooling 
in our experimental set-up (the shortest cooling time is 
about 30 minutes), we focus on the resulting fiuctuations 
(noise) in cr(t). After describing the samples and mea- 
surement techniques in more detail (Sec. [Ill), we present 
the analysis of the noise power spectra obtained using 
a AT protocol (Sec. IIIip . The results demonstrate the 
glassy arrest of the charge dynamics for < Tig, where 
the glass transition density rig is the same as that found 
in noise studies that employed An^ protocols. ^^^^ 

The power spectrum, being the Fourier transform of 
a temporal correlation function, only provides informa- 
tion about the second moment of the fiuctuating quan- 
tity. Additional information can be extracted from the 
full probability distribution of the fiuctuations. There- 
fore, in Sec. lIVl we compare the probability density func- 
tions (PDFs) obtained in the two protocols, Aug and 
AT. While in both cases the PDFs are non-Gaussian for 
Us < Tig, the change of the carrier density has a quali- 
tatively different and more dramatic effect than AT. In 
particular, the results suggest that Aug induces strong 
changes in the free energy landscape of the system, which 
should have important implications for theoretical mod- 
eling of the glassy dynamics in a 2DES. 

For a AT protocol, the power spectra and the second 
spectra^^^^ of the noise in the glassy regime have also 
been investigated for different values of the excitation 
voltage Vexc used in the measurements of a. In indium 
oxide, another realization of a Coulomb glass, the non- 
Gaussian noise was attributed to a non-linear mechanism 
caused by the high applied bias^^ As shown in Sec. |Vl 
in the 2DES in Si the large, non-Gaussian noise does not 
depend on the applied bias over more than two orders 
of magnitude variation of Kxc, providing additional evi- 
dence that it refiects the intrinsic glassiness of the 2DES. 
The main results are summarized and discussed further 
in Sec. |VD 



II. EXPERIMENT 

A. Samples and measurement techniques 

Measurements were carried out on rectangular n- 
channel (lOO)-Si MOSFETs with a relatively large 
amount of disorder. Similar to previous studies ^i^i"— 
the back-gate bias of —2 V was applied to maximize the 
4.2 K peak mobility to ~0.06 m^/Vs. The sample dimen- 
sions L X W (T-length, I^-width) were 1 x 90 fim^ for 
samples A and Al, and 2 x 50 /am^ for samples B and 



Bl. The devices were fabricated with poly- Si gates, self- 
aligned ion-implanted contacts, substrate doping Na ~ 
2 X 10^^ cm~^, oxide charge A^ox 1 x 10^^ cm~^, and 
oxide thickness dox = 50 nm. All four samples exhibited 
the same behavior, and some of them were, in fact, used 
in previous studies (e.g. A and B in Refs. flll - [l3 . Al in 
Ref. 9). 

The conductivity a was measured using a standard 
two-probe ac technique (typically at ~ 11 or 13 Hz) 
with an ITHACO 1211 current preamplifier, and either 
a PAR124A or a SR7265 lock-in amplifier. The con- 
tact resistances and the contact noise were determined 
to be negligible in these samples.^ Unless noted other- 
wise (see Sec. |V]), the excitation voltage ]4xc was kept 
constant and low enough (2 — 10 /iV) to ensure that the 
conduction was Ohmic. The low- frequency of ]4xc was 
limited by the low cut-off frequency of RC filters, which 
were used to reduce external electromagnetic noise, as 
well as by high resistance of the samples. A precision 
dc voltage standard (EDC MV116J) was used to apply 
the gate voltage Vg, which controls Ug'. ns(10^^cm~^) = 
4.31 (F^[V] — 6.3). Measurements of current (/) fiuctua- 
tions as low as 10~^^ A were performed both in a dilution 
refrigerator and in a ^He system (base T = 0.24 K). Rel- 
atively small fiuctuations of T, Vg, and Vq^c were ruled 
out as possible sources of the measured noise, since no 
correlation was found between them and the fiuctuations 
of /. The background noise was measured for all Ug and 
T by setting Vq^c to zero. It was always independent 
of frequency (i.e. white) and usually several orders of 
magnitude smaller than the sample noise. 



B. Temperature dependence of the time-averaged 
conductivity 

Figure [U shows the time-averaged (see Sec.lIIip conduc- 
tivity (a) as a function of T for different Ug for one of the 
samples. The density was varied at T ^ 10 — 12 K, where 
the 2DES is in a thermal equilibrium.^^ The critical den- 
sity Tie is determined from the fits to {a{ng,T)) on both 
metallic and insulating sides of the MIT, as discussed in 
detail elsewhere ^^li^ii^ii^ In particular, in the insulating 
regime, the transport is thermally activated, whereas in 
the intermediate, ric < rig < rig glassy phase, {a{ng,T)) 
acquires a very specific, non- Fermi liquid T^/^ correction. 
In all samples, ric = (4.5 ±0.4) x 10^^ cm~^, and the glass 
transition density determined from Arig protocols^iii"— 
was rig = (7.5 ±0.3) x 10^^ cm~^. The error bars account 
for sample-to-sample variations. 



III. CONDUCTANCE NOISE FROM A AT 
PROTOCOL 

In a AT protocol, each noise measurement starts by 
fixing the carrier density n^atT^lO — 12K. The sam- 
ple is then cooled to the measurement temperature. The 
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FIG. 1: (Color online) Sample B. Conductivity (a) vs. T for 
different ris, as shown. Us was varied at high T^^IO — 12K. 
Solid lines guide the eye. 



cooling times varied from 30 minutes to 10 hours, with 
no effect on the noise results. In addition, such slow 
cooling did not result in any visible relaxations, at least 
on our experimental time scales. Typical data are pre- 
sented in Fig.[2fa), which shows the relative fluctuations 
Aa{t) I (a) vs time for several Ug measured at T = 0.24 K. 
It is obvious that, by decreasing n^, the amplitude of the 
fluctuations increases dramatically. Even more striking 
is the change in the character of the noise: at low ng, fluc- 
tuations occur on many different time scales, including 
slow changes over periods of several hours; at high n^, on 
the other hand, the variance of the noise no longer varies 
with time. These features are more apparent in Fig.[2fb), 
where we choose to plot the same A(j{t) = {(j{t) — {a)) 
normalized by 6a, the square root of the variance of the 
fluctuations. The increase in the noise magnitude with 
decreasing ris is clearly accompanied by the onset of non- 
Gaussian behavior, similar to the noise results obtained 
from a An^ protocoL-iiSiii^ 

The normalized power spectra of the relative changes 
in conductivity Aa{t)/{(j) are found to obey the empiri- 
cal law Si (X (/ - frequency). In order to compare 
the noise magnitudes under different conditions, the spec- 
tra were averaged over octaves, and the resulting frac- 
tion of power Si{f = 1 mHz) is taken as the measure of 
noise. Figure [3l^ a) shows that Si{f = 1 mHz) increases 
by several orders of magnitude as is reduced below 
^ 7 X 10^^ cm~^, while a rises rapidly to a much higher 
value [Fig.[3l^b)] in the same range of n^. Both effects in- 
dicate that the electron dynamics suddenly and dramat- 
ically slows down for < rig ~ 7 x 10^^ cm~^, which is 
precisely the regime where non-Gaussian behavior of the 
noise is also observed (e.g. Fig. [2]). Therefore, the results 
of the noise measurements obtained using a AT protocol 
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FIG. 2: (Color online) Sample B. (a) Relative fluctuations 
Aa{t)/ (a) = {cr(t) — (cr)) / (a) vs time for several Us shown on 
the graph (top to bottom) at T = 0.24 K. The averaging of 
a{t) was done over a time interval shown on the plot. Traces 
are shifted for clarity. Us was changed at T?^ 10 — 12 K 
and the cooling time was about 10^ s. (b) In order to make 
the change in the character of the noise with Us more appar- 
ent, the same fluctuations Aa{t) = {cr{t) — (a)) are shown 
normalized by 6a, where 6a = ~ • 



(black triangles in Fig. [3j) also provide evidence for glassy 
freezing of the 2DES for < n^, with the value of Ug in 
agreement with that obtained from An^ protocols. 

Furthermore, the values of Si{f = 1 mHz,ns) and 
aijis) have been compared to those from a Aug pro- 
tocol (see, e.g.. Fig. |4]). In the latter, was changed 
in small steps (typically ~ 0.2 x 10^^ cm~^) at low T, 
where the equilibration times exceed the experimental 
time scales Figure [3] includes the results for Si and a 
obtained by changing at 1 K, 0.8 K (as in Ref. 0) and 
0.24 K. Both the noise magnitude and the exponent a 
demonstrate that AT and small Aug protocols produce 
the same results. 
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FIG. 3: (Color online) Open symbols: sample Al; solid sym- 
bols: sample B. Different symbols correspond to different 
T, as shown, where Us was changed before cooling down to 
0.24 K to measure noise, (a) The normalized noise power 
Si{f = 1 mHz) vs Us. (b) The exponent a vs Us. 



While both protocols, cooling with a fixed rig and a 
small change of Ug at low T, provide evidence for the 
onset of glassy dynamics from the fluctuations of cr(t), 
we note that there are no visible relaxations of a follow- 
ing either process. In both cases, though, the size of the 
perturbation {ks^T ~ 10 K or AEf ~ 1 — 2 K, respec- 
tively, where Ep is the Fermi energy) is much smaller 
than Ef (~ 20 — 60 K),^^ i.e. probably too small to lead 
to any observable relaxations. The relaxations were in- 
deed observed for much larger changes of ng,^^^^ such 
that AEf 30 - 150 K, i.e. for AEf > Ef ksT. 

Since AEf > ksT in all An^ protocols on a 2DES in 
Si MOSFETs, the perturbative treatments of the effect 
of the gate volt age, ^^^^^ which may be relevant for other 
materials, are not suitable for this system. In particu- 
lar, here the applied AVg (i.e. AEf) are expected to 
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FIG. 4: (Color online) Sample B. Conductivity a vs time at 
T = 0.24 K for several ris. Traces are shifted for clarity. The 
data were obtained (a) after changing Us at 10 K, followed 
by a slow decrease of T for about 12 hours (AT protocol), 
and (b) after a small ris change (Aris ~ 0.22 x lO^^cm"^) at 
low T = 0.24 K {Aus protocol). In another run, (c) Us was 
changed at 1 K (Aus protocol); after staying at 1 K for two 
hours, T was lowered to 0.24 K in one hour, and a{t) was 
then measured. 



trigger major, collective rearrangements of the electron 
configuration. As shown below (Sec. lIV A|) . a comparison 
of the full distribution of the fluctuations obtained from 
the two protocols demonstrates that even a small Aus 
{AEf Ef) has a qualitatively different effect than 
AT, even though the results for the noise power spectra 
(Fig.[3|) and ((j(ns,T)) are the same within experimental 
error. 
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IV. DISTRIBUTION FUNCTIONS OF THE 
FLUCTUATIONS 

In systems that are in equilibrium and away from criti- 
cality, the probabihty density functions (PDFs) of fluctu- 
ations in global quantities are Gaussian, according to the 
central limit theorem, which is based on the hypothesis 
that the system may be decomposed into many uncorre- 
lated elements. When this condition is not satisfied, e.g. 
due to the divergence of the correlation length at criti- 
cality or slow, nonexponential decay of correlation func- 
tions in glassy systems, non-Gaussian PDFs have been 
observed in a number of disparate systems. 



A. PDFs from An^ and AT protocols 

Figure [5] shows some typical PDFs of the conductance 
(G) fiuctuations for < Ug measured in a An^ protocol. 
The density n^, lying in the metallic glass regime {ric < 
Us < rig), was fixed at T ^ 0.8 K after a small change 
Aus = 0.43 X lO^^cm-^ {AEp - 3 K<C - 41 K). The 
sample was then cooled down to T = 0.13 K, and G{t) 
was recorded over sixteen sequential intervals of about 35 
minutes each. The PDFs, shown in Fig.[5fa) for three out 
of sixteen intervals, are clearly not Gaussian, but rather 
have a complicated, multi-peaked structure that changes 
with time. 

These results indicate that, on experimental time 
scales, the 2DES "visits" only a small number of states, 
but that new states also become available with time. 
This is also refiected in the correlated "wandering" of 
the PDFs with time [Fig. [5fb)]. For > n^, on the 
other hand, the PDFs become Gaussian (not shown) on 
much shorter time scales. 

The PDFs obtained from a AT protocol (Fig. [6l left 
column) exhibit the same general behavior: they are non- 
Gaussian for ris < rig and Gaussian for ris > rig. (Even 
though the non-Gaussian shape of the PDFs for ris < rig 
is more obvious on a semi- log scale, they are presented on 
a linear scale for easier comparison to the effects of An^. 
See also Sec. |Vl) Here, however, even the non- Gaussian 
PDFs are smooth, single-peaked functions, reminiscent 
of PDFs in a variety of systems displaying critical f^^i^i 
glassy^ or other out-of-equilibrium behavior (e.g. the 
Danube water level^^). In all these systems, the PDFs 
are skewed, resembling a zero-centered Gaussian, which 
describes pseudo-equilibrium fiuctuations, with one (ex- 
ponential) tail that is due to large, rare events. In that 
context, the totally irregular structure of the PDFs ob- 
tained after a small Arig (Fig. [5]) is striking. 

The above differences can be understood in the fol- 
lowing way. After cooling, the 2D charge glass settles 
in a deep energy valley, perhaps a global minimum of 
the FEL, corresponding to a given n^. As a function of 
time, the system will explore the sub valleys of the energy 
landscape, but it will in general remain in the "main", 
deep valley. Such pseudo-equilibrium fiuctuations will 




G-<G> (arb. units) 

FIG. 5: (Color online) Sample Al. Probability density func- 
tion (PDF) of the conductance (G) fluctuations for ric < 
ns(10^^cm"^) = 5.58 < at T = 0.13 K. ris was fixed at 
T 0.8 K after a smafl change An^ = 0.43 x 10^^ cm"^, cor- 
responding to AEf 3 Ef ~ 41 K. The measurement 
time was divided into sixteen intervals of about 35 minutes 
each. PDFs for intervals 1, 4, and 12 are shown in (a), demon- 
strating the non-Gaussian behavior of the noise, (b) The color 
map of the PDF of the fluctuations in C as a function of all 
sixteen sequential time intervals. The averaging (...) was per- 
formed over the total measurement time (^ 9 hours). 

produce PDFs of the type shown in the left column of 
Fig. [6l In the Arig protocol, in contrast, a change in the 
number of charge carriers reshuffles all energies because 
of the Coulomb interactions, so that the 2DES in general 
no longer finds itself in the minimum of the FEL. This 
will result in a more dramatic wandering of the system 
through the FEL, producing the PDFs of the type shown 
in Fig. [51 

B. PDFs from a waiting- time protocol 

Large An^ {AEp > Ep) lead to observable relaxations 
in cr(t). This effect has been used to investigate aging 
and memory in a 2DES in Si by employing the so-called 
waiting-time {tw) protocol, ^^^^ which consists of the fol- 
lowing. It starts with a AT protocol, i.e. cooling at 
a fixed density no and measuring noise (see left column 
in Fig. [6] for the corresponding PDFs). The density is 
then switched rapidly (within 1 s) to a different value 
ni, where it is kept for a time t^. During that time, the 
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FIG. 6: (Color online) Sample A. PDF vs [a{t) - (cr(t))]/(cro) 
at T = 0.24 K for several no, as shown, (ao) is the time- 
averaged conductivity corresponding to no at the measure- 
ment T, obtained from a AT protocol. The red curves are 
fits to a Gaussian distribution. Left column: PDFs of the 
noise after cooling from 10 K to 0.24 K for each given no 
(AT protocol). The cooling time was at least one hour long. 
In this protocol, (cr(t)) = (ao), i.e. there are no observable 
relaxations after cooling. Right column: PDFs of the noise 
measured after a subsequent change of Us from no to a much 
higher value m = 20.26 x 10^^ cm"^ during t^, = 1000 s. In 
this protocol, (cr(t)) describes the slowly relaxing background. 



system relaxes away from its initial (pseudo-) equilibrium 
state determined by no and towards a new equilibrium 
state determined by ni. Finally, the density is changed 
back (within 1 s) to no at t = 0, and the slowly evolving 
a{t) is measured. We note that here t = is defined 
as the time when the charge carrier density reattains its 
original value no- If tu, is shorter than the equilibration 
time, then a{t) will depend on t^^ i.e. on the measure- 
ment history (aging efect). It is also said that the system 
has a memory of the time it spent with ni. Here we focus 
on that situation, i.e. on noise during aging. 

Figure [6] presents a comparison of the PDFs obtained 
from AT (left column) and (right column) protocols 
for different no and the same ni = 20.26 x lO^^cm"^ and 
tw = 1000 s. For the highest no ^ n^, the PDFs are 
Gaussian in both protocols, with roughly the same stan- 
dard deviation. For lower no, however, there are some 
remarkable differences. First, the noise after a tempo- 



rary density change is much larger than the noise before 
the change, that is, after a AT protocol. At the lowest 
no, this difference amounts to more than one order of 
magnitude, as first noted in Ref. 13. Second, the change 
of density during results in complex, multi-peaked, to- 
tally random-looking PDFs, analogous to those obtained 
from An^ protocols (Sec. IIV Ap . These differences are 
even more striking considering that, for a given no, the 
PDFs in Fig. [6] were obtained from the measurements 
that were carried out under exactly the same experimen- 
tal conditions; the only difference was in the sample his- 
tory. 

These results imply the existence of a rugged free en- 
ergy landscape that is modified in a non-trivial way by all 
Aus of practical interest (corresponding to AEf > ksT). 
The rearrangements of the electron configuration and the 
reshuffling of the FEL caused by Coulomb interactions 
will clearly need to be taken into account in theoreti- 
cal models of the glassy dynamics of the 2DES in Si. 
This complicated problem is further aggravated by the 
proximity to the MIT, where substantial changes in the 
screening of the 2DES may be expected. ^^^^ 

The large noise that accompanies the aging process 
reflects the collective wandering of the 2DES through 
the FEL, as the system relaxes towards a global mini- 
mum that corresponds to the equilibrium state. In other 
out-of-equilibrium systems (e.g. polymers and colloidal 
glasses, conventional spin glasses^), the PDFs have 
been observed to evolve towards a Gaussian shape as the 
system ages. Here the evolution of the PDFs has been in- 
vestigated from the data obtained in several consecutive 
time intervals, during the relaxation of a{t) induced using 
a tw protocol, as described above. A typical behavior is 
presented in Fig. [71 which shows the PDFs for the same 
density no = 3.02 x lO^^cm"^ in the insulating regime 
(no < Tie) and for three different tw In each case, the 
PDFs for four consecutive 3-hour time intervals demon- 
strate that the noise indeed becomes more Gaussian with 
time, as the system slowly approaches equilibrium. 

Understanding of the age dependence of the PDFs is a 
non-trivial problem in general. Even in situations where 
PDFs feature a "simple" zero-centered Gaussian with an 
exponential tail, different models have been proposed. 
In addition, numerical simulations of aging models have 
been performed, including those based on energy land- 
scapes with hierarchically organized metastable states 
(see, e.g., Ref. 0). While it is unclear whether any 
of that work is relevant to our system, we note that re- 
sistance noise measurements on both spin glasses^^ and 
the Coulomb glass in a 2DES in Si^^ have provided sup- 
port for the hierarchical picture of glassiness.^^ A number 
of other, remarkable similarities in the phenomenological 
behavior of these two types of glassy systems have been 
also reported. iii — 
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FIG. 7: (Color online) Sample A. PDF vs. [a{t) - {a{t))]/{ao) 
at T = 1 K for no = 3.02 x lO^^cm"^, following a change 
of Us to m = 20.26 X lO^^cm"^ during (a) t^, = 10 s, (b) 
tw = 1000 s, and (c) tw = 31600 s. For each t"w 1 ^0 was 
initially set at 10 K. Different curves correspond to PDFs for 
four consecutive 3-hour time intervals, as shown. The PDFs 
become more Gaussian with time, as the system ages. 



V. DEPENDENCE OF CONDUCTANCE NOISE 
ON APPLIED ELECTRIC FIELD 



The emergence of the non-Gaussian conductance noise 
for Us < Tig^ i.e. in precisely the same regime where re- 
laxations exhibit the out-of-equilibrium behavior 
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weh as the detailed study and comparison of the noise ob- 
tained using different protocols (Sees. Hill and HV|). clearly 
demonstrate that it reflects the intrinsic glassiness of the 
2DES. For example, the qualitatively different behavior 
of the PDFs measured under the same conditions, but 



with different sample histories (Fig. [6j), rules out any 
spurious effects as possible sources of non-Gaussian noise 
from a protocol. Likewise, the agreement of the results 
for the noise power spectra obtained using AT and An^ 
protocols (Fig. [3]) implies that they have the same origin. 
Experiments on a disordered indium oxide, an electron 
glass deep in the insulating regime, have shown that,^^' 
in that system, the degree of non- Gaussian behavior may 
be affected by the applied electric bias Vq^c- In particu- 
lar, the non-Gaussian nature of the noise was found to be 
non-monotonic with lexc, vanishing at both low and high 
bias, and peaking at an intermediate value of Vq^c This 
behavior was attributed to a nonlinear effect inherent to 
variable-range hopping transport. While the indium ox- 
ide samples were indeed very deep in the hopping regime, 
this is in general not true of the experiments on a 2DES, 
where glassiness sets in already on the metallic side of the 
MIT {ug > Tie). Nevertheless, it is interesting to investi- 
gate whether similar effects might be observed also in the 
glassy regime of the 2DES, in particular in the range of 
parameters where non-Gaussian conductance noise has 
been studied (Sees, [ml [iVland Refs. 9^10^5). 

In order to explore the effects of ]4xc, conductance 
noise Aa{t)/{a) was measured using a AT protocol. Af- 
ter changing Ug at a high T (here at ^ 20 K), the 2DES 
was cooled to a desired T, and a{t) measured for constant 
values of V^xc ranging from 0.5 /iV to 1 mV. Typical data 
are presented in Fig. [51 which shows Act/ (<j) vs time for 
Us = 3.24 X 10^^ cm~^ in the insulating regime (n^ < Uc) 
and for several Vexc at two different T. The noise is obvi- 
ously non-Gaussian and large, of the order of 100% at low 
T. The amplitude of the relative fluctuations decreases 
with increasing T, as established in earlier studies, ^^^^ as 
well as for high enough excitations. 

It is indeed obvious from the corresponding normalized 
power spectra Si{f) = C/f^ (Fig.[9j) that the noise does 
not depend on the applied excitation at low enough lexc, 
but the power spectra Sj become suppressed as Vq^c is 
increased beyond a certain value. A detailed comparison 
of the dependence of the exponent a, the parameter C, 
and the time-averaged conductivity (a) on V^xc is pre- 
sented in Fig. [TOl It is apparent that {a) remains Ohmic 
over a wide range of Vq^c [Fig. [TOlf a)], up to an electric 
field - 50 V/m (Kxc ^ 100 /iV) at T - 0.1 K (where 
eVexc ^ kBT). In the same range of Kxc, both the ex- 
ponent a ^ 1.6 [Fig. [TOlf b)] and the noise magnitude, i.e. 
the parameter C [Fig. [TOlf c)]. also do not depend on V^xc 
within experimental error. The latter result is expected 
if the applied voltage is sufficiently small to be in the 
Ohmic regime and if the noise is in fact due to conduc- 
tance (or resistance) fluctuations. At higher 14xc, where 
(cr) deviates from the Ohmic behavior, a does not seem 
to change with 14xc or T within the scatter of data, but 
the parameter C clearly decreases rapidly with increas- 
ing T4xc5 similar to the effect of raising T. Therefore, the 
results strongly suggest that the suppression of the noise 
by high V^xc is most likely caused by simple heating. 

Similar to previous studies of noise on the 2DES in 
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FIG. 8: (Color online) Sample Bl. Relative fluctuations 
of conductivity (cr(t) — (a)) /{a) vs time for Us = 3.24 x 
10^^ cm"^ < Tic and different Kxc at a) T = 0.100 K and 
h) T — 0.300 K. Traces are shifted for clarity. The data were 
obtained after changing ris at T 20 K, cooling down to the 
desired T, and then measuring a{t) for different values of Kxc- 



S i-"J"^ and those in indium oxide^^, the effect of Ve^c on 
the non- Gaussian character of the noise has been inves- 
tigated by analyzing the second spectrum 5*2 (/2,/), a 
fourth-order noise statistic. ^^^^ In particular, 5'2(/2,/) 
represents the power spectrum of the fluctuations of 
Si if) with time. For Gaussian noise, 5*2 (/2) is white (in- 
dependent of /2). On the other hand, for non-Gaussian 
noise, 5*2 ex l//2~^ with the exponent (1— /3) ^ 0. There- 
fore, the value of the exponent (1 — /3) is a convenient 
measure of the non- Gaussian character of the noise. S2 
was analyzed by using digital filtering^^'^^ in a given fre- 
quency range / = {fh^fH) (usually fn = 2/l). Some 
examples of the normalized 6*2 (/2), with the Gaussian 
background subtracted, are presented in Fig. fTTT a) for 
the data in Fig. [8l^a). The second spectrum is clearly 
non- white, in agreement with earlier work.^^ Here we 
find that the exponent (1 — /3) ^ 0.8 does not depend 
on 14xc [Fig. fTTlb)] in the regime studied. Therefore, 
non- Gaussian noise is observed over a wide range of elec- 
tric fields in the regime of linear or Ohmic response. 

Detailed studies of non-Gaussian noise in the 2DES in 
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FIG. 9: (Color online) Sample Bl. The octave- averaged noise 
power spectra Si — C / \s f for several V^xc at a) T = 
0.100 K and b) T = 0.300 K. SoUd lines are linear least- 
squares fits with the slopes equal to a. 



Si (Sees. imillYland Refs. HUolIill) have been performed 
with the applied electric fields ^ 5 V/m, i.e. an order of 
magnitude lower than the field where the non- Ohmic be- 
havior sets in (see Fig. [TQ|). This is in contrast to the ex- 
periment on indium oxide, where only the non-Ohmic 
regime was explored, and where the non-monotonic de- 
pendence of the non-Gaussian behavior on Vq^c was re- 
ported. A study of such a deep non-hnear regime, how- 
ever, is beyond the scope of this work, so the existence of 
a similar effect in our system for much higher 14xc can- 
not be ruled out. Our study does demonstrate, on the 
other hand, that the non-Gaussian noise in a 2DES is 
not due to the same non-linear mechanism caused by the 
high bias. Instead, it provides further evidence that the 
correlated noise reflects the intrinsic out-of-equilibrium 
dynamics of the 2DES in Si. 



VI. CONCLUSIONS 

We have presented a study of the conductance noise in 
a 2DES in Si at low T, measured by using three different 
experimental protocols. In the so-called AT protocol, the 
carrier density is fixed at a high T where the system is 
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FIG. 10: (Color online) Sample Bl. (a) (a), (b) exponent a, 
and (c) parameter C (in Si = C//") vs Ve^c at T = 0.100 K 
(red squares) and T = 0.300 K (blue dots) for ris = 3.24 x 
10^^ cm~^ < ric- In (a) and (c), the arrows show the onset of 
the non-Ohmic regime, and dotted lines guide the eye. In (b), 
dotted lines represent the average values of a in the Ohmic 
regime: a = 1.5 db 0.2 at T = 0.100 K and a = 1.7 ± 0.1 at 
T 0.300 K. 



in an equilibrium state, and the noise is measured after 
cooling with a fixed n^. The results provide evidence 
that, similar to other types of glasses, the 2DES falls out 
of equilibrium as T is reduced. This glassy freezing is 
observed for all densities < n^, as T ^ 0. 

In a Aus protocol, Ug is changed by a small amount 
at such low T that the system is unable to equilibrate 
on experimental time scales. While the perturbations 
introduced by both AT and An^ protocols are too small 
to lead to observable relaxations of a (i.e. kB^T^ AEf ^ 
Ef)^ they produce the same results for the power spectra 
of the non-Gaussian noise observed for ng < n^, where 
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FIG. 11: (Color online) Sample Bl, T 



0.100 K, Us 



10^^ cm ^. (a) Second spectral density 6*2 oc I//2 



3.24 X 

vs /2 

for several values of V^xc, as shown. 6*2 was measured in the 
frequency octave / = (/l,2/l), where /l = 16 mHz. The 
solid lines are linear least-squares fits with the slopes equal to 
1 — f5. (b) Exponent (1 — /3) vs Ksxc for different values of /l, 
as shown. Dotted line guides the eye. 



the glass transition density Ug is found to be the same 
in both cases. The analysis of the probability density 
functions (PDFs) of the fluctuations, however, reveals 
that Aus has a strikingly different effect on the system 
than AT (both AEF.ksAT > ksT). In particular, the 
results strongly suggest that the density change reshuffles 
all energies, because of the Coulomb interactions, thus 
modifying the free energy landscape of the 2DES. For 
this reason, theoretical modeling of the glassy dynamics 
in this system might be considerably more difficult than 
in some other types of glassy materials. 

The effect of the density change was studied further by 
using a waiting-time (tw) protocol, in which Us is changed 
only temporarily during t^^ but AEf > Ef- As a result 
of such a large perturbation, the system exhibits visible 
relaxations of a, i.e. aging.—"— Interestingly, the non- 
Gaussian PDFs exhibit both history dependence and an 
evolution towards a Gaussian shape as the system ages 
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and slowly approaches equilibrium, similar to the behav- 
ior of a great variety of out-of-equilibrium systems. 

The power spectra and the higher order noise statistics 
of the non- Gaussian noise have been also investigated 
over a wide range of the applied bias voltage. The data 
demonstrate that the non-Gaussian noise is observed in 
the regime of linear response, i.e. that it is not caused by 
the application of a high bias, but rather that it reflects 
the intrinsic out-of-equilibrium behavior of the 2DES. 

In summary, we have established several new char- 
acteristics that the 2DES in Si has in common with a 
large class of both 2D and 3D out-of-equilibrium systems. 
Thus this work provides additional strong evidence that 
many such universal features are robust manifestations 
of glassiness, regardless of the dimensionality of the sys- 



tem. In addition, we have revealed some effects that 
are unique to Coulomb glasses. Therefore, our findings 
should be helpful in the understanding of the complex 
behavior near the MIT in a variety of strongly correlated 
materials. 



Acknowledgments 

This work was supported by NSF DMR-0905843, the 
National High Magnetic Field Laboratory through NSF 
Cooperative Agreement DMR-0654118, and the State of 
Florida. 



* Electronic address: lin@ magnet. fsu.edu] 

^ Electronic address: shi@magnet.fsu.edu] 

* Electronic address: j aroszy@ magnet . f s u^ edu | 

^ Electronic address: dragana@ magnet . f su . edu ' 

^ V. Dobrosavljevic, in Conductor- Insulator Quantum Phase 

Transitions, edited by V. Dobrosavljevic, N. Trivedi, and 

J.M. Valles Jr. (Oxford University Press, 2012, ISBN 

9780199592593). 
^ E. Miranda and V. Dobrosavljevic, Rep. Prog. Phys. 68, 

2337 (2005). 

^ E. Miranda and V. Dobrosavljevic, in Conductor- Insulator 
Quantum Phase Transitions, edited by V. Dobrosavljevic, 
N. Trivedi, and J.M. Valles Jr. (Oxford University Press, 
2012, ISBN 9780199592593). 

^ L. C. E. Struik, Physical Aging in Amorphous Polymers 
and Other Materials (Elsevier, Amsterdam, 1978). 

^ See Complex Behavior of Glassy Systems, edited by M. 
Rubi and C. Perez- Vicente (Springer, Berlin, 1997). 

^ For a review, see Dynamic Heterogeneities in Classes, Col- 
loids and Cranular Media, edited by L. Berthier, G. Biroli, 
J.-P. Bouchaud, L. Cipelletti, and W. van Saarloos (Oxford 
University Press, 2011). 

^ For a review, see A. Amir, Y. Oreg, and Y. Imry, Annu. 
Rev. Condens. Matter Phys. 2, 235 (2011). 

® D. Popovic, in Conductor- Insulator Quantum Phase Tran- 
sitions, edited by V. Dobrosavljevic, N. Trivedi, and 
J.M. Valles Jr. (Oxford University Press, 2012, ISBN 
9780199592593). 

^ S. Bogdanovich and D. Popovic, Phys. Rev. Lett. 88, 
236401 (2002). 

J. Jaroszyiiski, D. Popovic, and T. M. Klapwijk, Phys. 
Rev. Lett. 89, 276401 (2002). 

J. Jaroszyiiski and D. Popovic, Phys. Rev. Lett. 96, 037403 

(2006) . 

J. Jaroszynski and D. Popovic, Phys. Rev. Lett. 99, 046405 

(2007) . 

J. Jaroszynski and D. Popovic, Phys. Rev. Lett. 99, 216401 
(2007). 

J. Jaroszynski and D. Popovic, Physica B 404, 466 (2009). 
J. Jaroszynski, D. Popovic, and T. M. Klapwijk, Phys. 
Rev. Lett. 92, 226403 (2004). 

V. Dobrosavljevic, D. Tanaskovic, and A. A. Pastor, Phys. 



Rev. Lett. 90, 016402 (2003). 

G. F. Rodriguez, G. G. Kenning, and R. Orbach, Phys. 
Rev. Lett. 91, 037203 (2003). 

M. B. Weissman, Rev. Mod. Phys. 60, 537 (1988). 

M. B. Weissman, N. E. Israeloff, and G. B. Alers, J. Magn. 

Magn. Mater. 114, 87 (1992). 

M. B. Weissman, Rev. Mod. Phys. 65, 829 (1993). 
V. Orlyanchik, V. I. Kozub, and Z. Ovadyahu, Phys. Rev. 
B 74, 235206 (2006). 
22 ^^[K] = 7.31n,[10^^ cm-2] for electrons in Si 
MOSFETs.22 

T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 
437 (1982). 

E. Lebanon and M. Miiller, Phys. Rev. B 72, 174202 
(2005). 

M. Miiller and E. Lebanon, J. Phys. IV (France) 131, 167 
(2005). 

S. T. Bramwell, P. C. W. Holdsworth, and J.-F. Pinton, 
Nature 396, 552 (1998). 

S. Joubaud, A. Petrosyan, S. Ciliberto, and N. B. Garnier, 
Phys. Rev. Lett. 100, 180601 (2008). 

S. T. Bramwell, T. Fennell, P. C. W. Holdsworth, and B. 
Portelh, Europhys. Lett. 57, 310 (2002). 
A. A. Pastor and V. Dobrosavljevic, Phys. Rev. Lett. 83, 
4642 (1999). 

M. Miiller and L. B. loffe, Phys. Rev. Lett. 93, 256403 
(2004). 

^'^ S. Pankov and V. Dobrosavljevic, Phys. Rev. Lett. 94, 
046402 (2005). 

L. Buisson, M. Ciccotti, L. Bellon, and S. Ciliberto, in 
Fluctuations and Noise in Materials, edited by D. Popovic, 
M. B. Weissman, and Z. A. Racz, Proceedings of SPIE 
5469 (SPIE, Bellingham, WA, 2004), pp. 150-163. 
P. Sibani, G. F. Rodriguez, and G. G. Kenning, Phys. Rev. 
B 74, 224407 (2006). 

A. Fischer, K. H. Hoffmann, and P. Sibani, Phys. Rev. E 
77, 041120 (2008). 

K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 
(1986). 

2^ G. T. Seidler and S. A. Solin, Phys. Rev. B 53, 9753 (1996). 
K. M. Abkemeier, Phys. Rev. B 55, 7005 (1997). 



